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ABSTRACT 

We have performed new wide-field photometry of the young open cluster NGC 6231 to study 
the shape of the initial mass function (IMF) and mass segregation. We also investigated the 
reddening law toward NGC 6231 from optical to mid-infrared color excess ratios, and found that 
the total-to-selective extinction ratio is Ry = 3.2, which is very close to the normal value. But 
many early-type stars in the cluster center show large color excess ratios. We derived the surface 
density profiles of four member groups, and found that they reach the surface density of field 
stars at about 10', regardless of stellar mass. 

The IMF of NGC 6231 is derived for the mass range 0.8 - 45 M Q . The slope of the IMF 
of NGC 6231 (r = —1.1 ± 0.1) is slightly shallower than the canonical value, but the difference 
is marginal. In addition, the mass function varies systematically, and is a strong function of 
radius - it is is very shallow at the center, and very steep at the outer ring suggesting the cluster 
is mass segregated. We confirm the mass segregation for the massive stars (m > 8 M Q ) by a 
minimum spanning tree analysis. Using a Monte Carlo method, we estimate the total mass of 
NGC 6231 to be about 2.6(±0.6) x 10 3 M . We constrain the age of NGC 6231 by comparison 
with evolutionary isochrones. The age of the low-mass stars ranges from 1 to 7 Myr with a slight 
peak at 3 Myr. However the age of the high mass stars depends on the adopted models and is 3.5 
± 0.5 Myr from the non- or moderately-rotating models of Brott et al. as well as the non-rotating 
models of Ekstrom et al. But the age is 4.0 - 7.0 Myr if the rotating models of Ekstrom et al. are 
adopted. This latter age is in excellent agreement with the time scale of ejection of the high mass 
runaway star HD 153919 from NGC 6231, albeit the younger age cannot be entirely excluded. 

Subject headings: stars: formation - stars: pre-main sequence - open clusters and associations: individual 
(NGC 6231) 
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1. INTRODUCTION 

The young open cluster NGC 6231 is consid- 
ered as the core cluster of the Sco OBI associa- 



tion (|Perrv et al.l Il990f) and is one of the bright- 
est open clusters in the southern hemisphere. Al- 
though there is a small amount of differential red- 
dening (|Sung et al.lll998HRaboud et al.lll997t ). the 
reddening of the inner 10' region is relatively con- 
stant and low. NGC 6231 is a rare case of an 
exposed young open cluster. In addition, there 
are about 150 O, B, and A main sequence (MS) or 
post-MS stars within 10' from the center of NGC 
6231, and therefore NGC 6231 is a good target 
for t he study of dynamical evolut ion of open clus- 
ters (jRaboud fe Mermilliodlll998l ). In most cases 
there are only a handful of massive O-type stars 
in nearby open clusters (the Orion Nebula Clus- 
ter or NGC 2264) and as a result, the initial mass 
function (IMF) of the massive part of these clus- 
ters inevitably suffers from low-number statistics. 
NGC 6231 is richer in massive stars than other 
nearby young open clusters, and is relatively closer 
and more transparent t han young compact c lus- 
ters such as NGC 3603 (|Sung fc Bessel]ll2004h . In 
this respect, NGC 6231 is a good target for the 
study of the stellar IMF, especially in the massive 
part. 

One of hot issues in cluster research is the ori- 
gin of mass segregation in open clusters. The 
fa ilure of the classica l appr oach was well tested 



Moeckel fc Bonnelll (|2009l ). and they concluded 



that several young open clusters, including the 
Orion Nebula Cluster and NGC 6231, are too 
young to be dynamically relaxed. They also con- 
cluded that strong primordial mass segregation is 
unlikely to be dynamically erased over the first few 
Myr of cluster evolution. The breakthrough in the 
so-called "primordial mas s segregation" in youn g 
open clusters was made bv lMcMillan et al.1 (120071 ). 
They showed that mergers of small clumps, that 
are initially mass segregated, lead to larger sys- 
tems th at inherit the p rogenitor clu mps' segre- 
gation. lAllison et al~l (|2009bl l2010h performed 
numerical simulations of cool, clumpy clusters, 
and showed that subvirial, fractal clusters are of- 
ten dynamically mass segregated on timescales far 
shorter than that expected from simple dynami- 



cal models. Allison et al. 1 (l2010l) showed that the 
formation of Trapezium-like massive multiples is 
highly probable in such a cool, clumpy cluster at 
a very early phase. They also showed the possi- 
bility of ejection of massive stars due to the dy- 
namical interaction between massive stars during 
the c ore collapse p h ase of the cluster. In addi- 
tion, Allison et al. ( 2009al ) devised an objective 
method to calculate the degree of mass segrega- 
tion based on the minimum spanning tree method. 
One of the important aims of this work is to study 
the extent of mass segregation in NGC 6231. 

Sung et all (|l998h (hereafter SBL98) performed 



UBVRI & Ha photometry and found several pre- 
main sequence (PMS) stars and PMS candidates. 
But the number of PMS stars with Ha emis- 
sion was far smaller than the number of low-mass 
PMS stars expected from the number of mas- 
sive MS and post-MS stars in the cluster. They 
suggested possible explanations for the deficit of 
PMS stars with Ha emission, such as past su- 
pernova explosions or strong stella r winds from 
massive stars in the cluster. Later Baume et al 



(1999) undertook a statistical approach and con- 
firmed the existence of a large number of PMS 
stars in the cluster. From analysis of the 180-ks 



XMM-Newton campaign of NGC 6231. lSana et al 
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( 20061 ) found 610 X-ray sources in the field of view 
and confirmed that many of these X-ray sources 
are bona-fide low -mass PM S stars in NGC 6231 
(jSana et al.ll2007l ). Recently iDamiani et all (|2009h 
detected more than 1600 X-ray point sources from 
the 120-ks Chandra ACIS o bservation, and con- 
firmed the ISana et al. ( 2007 ) results. In addition, 
they analyzed the X-ray spectra of about 500 X- 
ray sources without optical or near-IR counter- 
parts, and found unusually hard spectra for these 
sources. These objects could be low-mass PMS 
stars (V — I > 2.0) showing high coron al temper- 
atures (see Fig. 1 of ISung et al. ( 2008a )). 

Another interesting issue for NGC 6231 con- 
cerns past supernovae explosions in the cluster. 
The massive star content, the IMF, and the age 
of NGC 6231 strongly suggest the l ikelih ood of 
at least one supernova. iBessell et al. I ([l974h stud- 
ied the 06.5Iaf star HD 153919, the optical com- 
panion of the X-ray binary, 2U1700-37, and sug- 
gested from the reddening and brightness of the 
star, that it may be at th e same distance as Sco 
OBI. Later lAnkay et al.1 (|200ll ) suspected that 
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Sco OBI may be the origin of this runaway high- 
mass X-ray binary using the high-quality proper 
motion data from Hipparcos. From the optical po- 
larization observation of 35 stars in the cluster, 
Feinstein et al. I (l2003h found a curious semicircu- 
lar pattern of polarization angles of some stars in 
the center of NGC 6231. They interpreted it as 
a fingerprint of a past supernova explosion. But 
they seem not to have realized the existence of the 
runaway hi gh-mass X-ray bin ary HD 153919. Al- 
ternatively, |Bl^ck_et_alJ ( 2008 ) proposed the semi- 
circular pattern of polarization angles as the rem- 
nants of the historical nova in 1437 recorded in 
the old Korean Veritable Records of King Sejong 
Sejong Sillok. 

iRaboud (1996) studied the binary frequency of 
B type stars in NGC 6231 an d derived a mini- 
mum fraction of 52%. Recently ISana et al. ( 20081 ) 
studied the binary fraction of O type stars in 
the cluster and obtained a minimum fraction of 
0.63. The fraction of binaries, as well as the 
mass ratio in binary systems , could be very im - 
portant subjects of study (see Sana et al. ( 20121 )). 
because we still have insufficient information on 
the mass function o f secon daries in binary systems 
see IWeidner et ail (j2009l) for its influence on the 
mass function). In this respect NGC 6231 could 
be a good target because it is relatively closer and 
therefore it is easy to obtain high quality spectra 
with medium-size telescopes. 

New wide-field UBVI and Ha photometry of 
NGC 6231 is presented and compared in section 
2. Based on the new photometric data, several 
photometric diagrams are presented in section 3. 
The reddening law toward NGC 6231, reddening, 
distance, and surface density profiles are also de- 
rived in section 3. The Hertzsprung- Russell dia- 
gram (HRD) of NGC 6231 is constructed in section 
4. The IMF and age of NGC 6231 is derived from 
the HRD by comparing the observed data and stel- 
lar evolution models in this section. The mass seg- 
regation in NGC 6231 is investigated in the same 
section. Discussion on the runaway 06.5Iaf star 
HD 153919 and the origin of mass segregation is 
presented in section 5. The summary is in section 
6. 



2. OBSERVATIONS AND DATA RE- 
DUCTION 

UBVI and Ha CCD photometry of NGC 6231 
was performed on 2000 June 22 (NW & SW re- 
gions) and 25 (NE & SE regions) at Siding Spring 
Observatory with the 40 inch (1 m) telescope (f/8) 
and a thinned SITe 2048 x 2048 CCD (24 pix- 
els). The filters used were the same as those used 



for t he observation of NGC 2244 (jPark fc Sung 
20021 ). Two sets of exposure times were used in 
the observations - long: 60 s in /, 150 s in V, 300 
s in B, 600 s in U and Ha - and short: 2 s in / 
and V, 3 s in B, 30 s in U and Ha. The seeing 
was good on 2000 June 25 (about l."4 in V long- 
exposed images), but was bad on 2000 June 22 
(2. "2 - 3."2). The seeing strongly affects the pho- 
tometry of faint stars - completeness and limiting 
magnitude. The limiting magnitude in V on 2000, 
June 22 (Vn m u ~ 20.5 mag) is more than 0.5 mag 
brighter than that on 2000, June 25 {V Umtt > 21.0 
mag). 

All the preprocessing, such as the overscan 
correction, bias subtraction, and flat fielding, 
was done using the IRAF/CCDRED package. 
Instrumental magnitudes were obtained using 
IRAF/DAOPHOT via point spread function fit- 
ting for the target images and via simple aperture 
photometry for standard stars. All the instrumen- 
tal magnitudes were transformed to the standard 
UBVI system using SAAO E-region standard 
stars, E5 (E5-R, -O, -T), E6 (E6- U, 61, -Y, -W) 
and E 7 (E7-W, 37, -c, -X) (see ISung fc Besselll 
(|2000l ) for details). We observed E5 (and E7 on 
2000, June, 25) region(s) twice at high (airmass 

1.1) and low altitude (airmass ss 2.0) to deter- 
mine the atmospheric extinction coefficients and 
the time variation coefficients accurately. The 
primary extinction coefficients, time variation co- 
efficients of the photometric zero points (unit: 
magnitude per hou r - see equations (1) - (4) in 



Sung et al 



(20083)), 

and photometric zero points 
at midnight were derived from standard star pho- 
tometry, and listed in Table HJ On 2000 June 22, 
the time variation was a function of wavelength, 
largest in U but zero in /. On 2000 June 25, the 
time variation in U BV had similar values, but in 
I it had the opposite sign. The secondary extinc- 
tion coefficients for U an d B were adopt e d from 
the mean values used in Sung fc Besselll ( 2000l ). 
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Although there were no standard magnitudes in 
Ha, we transformed Ha to R ( i.e. pseudo-i?) so 



as to calculate th e Ha index ( Sung et al. 200(A 



Park fc Sund 120021) 



A total of 27,325 stars were measured in about 
a 40' x 40' area centered on NGC 6231, and are 
presented in Table [2] For 12 stars in Table [2] that 
were saturated in our images we have substituted 
data from SIMBAeH, WEBD^I (Open cluster 
database), or SBL98. Data for some stars missed 
in the new photometry were taken from SBL98. 
These data are marked with an "*" after the iden- 
tification number in the first column. We also 
listed 2MASS identification of th e objects in the 



2MA SS point source catalogue (jSkrutskie et al 



20061 ) in the 16th col umn. The optica l counter 



parts of X-ray sources ( Sana et al. 20061 ) are iden 
tilled, and marked as "X" (the closest or the only 
objects within the search radius of 3. "5) or "x" 
(star within 3."0, but second closest star). We 
also identified Ha emission stars (see section 3.1), 
and these are marked as "H" (strong Ha emission 
stars) or "h" (weak Ha emission candidates) . The 
membership information is in the 17th column. A 
total of 522 X-ray emission stars ("X" in column 
17) and 70 X-ray emission candidates ("x") are 



identified among 610 X-ray sources (jSana et al 
120061 ). Among them, 11 stars are X-ray emission 
stars with Ha emission ("+"), and 7 stars are X- 
ray emission stars with weak Ha emission ( "-" ) . 
Other information such as duplicity, variability, ID 
in SBL98, spectral type, and other identifications 
are listed in subsequent columns. The finder chart 
based on the new photometry is shown in Figure 

m 

The UBVI magnitude and colors from this 
photometry were compared with the previous pho- 
toelectric and CCD photometry in Table Al- 
though the agreement in zero points was quite 
good, the scatter was somewhat large because 
most photoelectric photometric studies were lim- 
ited to bright stars and many of these bright stars 
are variables or optical doubles. We compared our 
new data with SBL98 along with the relevant col- 
ors in Figure [2] The major difference between 
the two sets of data was in the standard stars 
used. Our new photometry, relative to SBL98, was 
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slightly brighter in V and slightly bluer in B — V 
for redder stars. Our new V — I was slightly red- 
der, but with no evident color dependency. For 
U — B, although the average difference is close to 
zero, there was a strong color-dependency. Such 
a difference is caused by the non-linear correc- 
tion term in t he transformation to t he Landolt's 
U system (see Sung k. Bessell ( 2000l ) for details) 
which is related to the size of the Balmer jump. 
In this respect (see the bottom panel of Figure [5]) 
we can understand the difference between our new 
U — B and those from photoelectric photometry, 
i.e. most stars observed with photoelectric pho- 
tometry are brighter and bluer, and therefore our 
new U — B is slightly bluer for blue stars. 

3. PHOTOMETRIC DIAGRAMS AND 
SURFACE DENSITY PROFILES 

3.1. Color-Color Diagrams and Reddening 

The color-color diagrams of the stars in the ob- 
served region are shown in Figure [3l X-ray emis- 
sion stars and Ha emission stars are marked by 
different symbols. The intrinsic and reddened MS 
relations (using the adopted reddening law for the 
(B — V, V — I) diagram) are also drawn in the 
figures. 

The (U — B, B — V) diagram is the basic di- 
agram for estimating the reddening to early-type 
stars. The E(B — V) reddening of early-type stars 
in NGC 6231 is between 0.45 and 0.60 mag, but 
some stars have somewhat larger E(B — V). The 
mean value of E(B — V) is 0.47 mag, which is 
the same as the value obtained by SBL98. The 
low-mass PMS stars with X-ray and/or Ha emis- 
sion seem to have nearly the same reddening as 
the massive early-type stars. But most of the 
normal late-type stars are slightly less reddened, 
and are presumably foreground field stars in the 
Sagittarius ar m. Unlike the low -mass PMS stars 
in NGC 2264 (jSung et al.lll997l) . only a few stars 
show ultraviolet (UV) excesses. This means that 
the accretion activity of PMS stars in NGC 6231 
has nearly stopped. We determined the reddening 
E{B — V) of all early-type stars (We assumed the 
standard color excess ratio of 0.72 between B — V 
and U — B), and calculated the spatial variation 
of reddening as shown in Figure |4j The dots in 
the figure represent the early type stars used in 
the reddening determination. As already noted in 
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SBL98, the reddening is very high in the south, 
relatively low in the north, and shows a local min- 
imum near the center. 

The (B — V, V — I) diagram is shown in the 
middle panel of Figure |3l The reddened MS rela- 
tion does not pass through the mean distribution 
of early-type stars. This is because the direction 
of the reddening vector in the diagram is affected 
by the reddening law and the adopted reddening 
law does not fit the early-type stars in the cluster 
center (see section 3.2 for details). 

The (Ha, V — I) diagram is presented for 
V < 18 mag in the lower panel of Figure [3J The 
low-mass PMS stars with Ha emission can be se- 
lected from the diagram ([Sung et al.l Il9971 ). As 
already noticed in SBL98, there are not many 
Ha emission PMS stars in NGC 6231. Most X- 
ray emission stars do not show excess emission in 
Ha. We designated stars above the border line 
(|Park fc Sunell2002h as Ha emission stars. In ad- 
dition, the Ha emission stars selected in SBL98 or 
from other studies are also designated as Ha emis- 
sion stars. A total of 31 Ha emission stars and 
34 Ha emission candidates were selected. Among 
them 18 stars are also X-ray emission stars with 
membership "+" or "-" . The star ID 357 (= HD 
326324 = Hen 3-1269) does not show any apprecia- 
ble emission in our Ha photometry, but t he star 
was c lassified as an Ha emission star by iHenize 
(ll976Tl . The stars ID 294 (= Wray 15-1546. (fWrav 
Il966h ) and ID 6017 (= Hen 3-1281) show strong 
Ha emission from our Ha photometry. ID 2132 
(=SBL 125) and ID 4359 (=SBL 667 = Se 166) 
are two massive Ha emission stars near the MS 
band within the cluster radius. 

3.2. Reddening Law 

The interstellar reddening law is one of the fun- 
damental parameters used in determining the dis- 
tance to astronomical objects, and is known to be 
different from sight line to s ight line in the Galaxy 
(|Fitzpatrick fe Massal [20091 ) . They confessed that 
there is no universal near-infrared (near-IR) ex- 
tinction law. We derived similar relations be- 
tween Ry and color excess ratio using their equa- 
tion (5) (the effect ive wavelength of Spitzer IRAC 
(|Fazio et al.ll2004l) bands from IRAC Instrument 
Handbook^) and power a from their Table 4. The 



results are 
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0.361, r = 0.998, (1) 



- 0.248, r = 0.999, (2) 



0.168, r = 1.000, (3) 



0.101, r = 1.000. (4) 



We also derived similar relations for the near- 
IR 2MASS bands, but the data points showed 
a large scatt er. We decided t o use the relations 
presented in iGuetter fc Vrbal (|l989t ) for near-IR 
JHK S bands. 

SBL98 obtained a slightly higher value of 
R v for the stars in NGC 6231 (R v = 3.3) us- 
ing the color excess in V — I. In this study 
we try to determine the reddening law from 
the optical I band to the mid-infrared (mid- 
IR) Spitzer IRAC bands. Firstly one of the au- 
thors (M.S.B.) calculated the synthetic colors of 
early-type stars in the near-IR 2MASS bands as 
well as in the mid-IR Spitzer IRAC bands. In 
addition, we downloaded the Spitzer IRAC im- 
ages from the Spitzer archivqj, and measured 
the magnitudes of O and B type stars in the 
same way as in ISung et all (|2009h . NGC 6231 
was observed as a part of the GLIMPSE (Galac- 
tic Legacy Infrared Mid-Plane Survey Extraor- 
dinaire) survey, but only part of the cluster was 
observed. Our photometric data are compared 
with those in the GLIMPSE Point Source Archive 
( htt p : //ir sa . ipac . caltech . edu/data/SP ITZER/GL IMPSE/[ ) , 
and the results are A[3.6] = +0.003 ± 0.132 
(n = 14), A[4.5] = -0.003 ± 0.071 (n = 21), 
A[5.8] = +0.038 ±0.115 (n = 14), and A[8.0] = 
+0.032 + 0.063 (n = 21). 

The color excess relations between B — V and 
other colors are presented in Figure [5] Although 
many early-type stars with E(B — V) = 0.4 - 
0.6 show a large scatter, if we use the stars with 
E(B — V) > 0.6, the color excess ratios are very 
close to the normal value. Using the relations be- 
tween Rv and color excess ratios, the total to se- 
lective extinction ratio Rv of the stars in NGC 
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6231 is 3.22 ± 0.04 from 8 optical to mid-IR col- 
ors. This value implies that the reddening law 
toward NGC 6231 is very close to normal. We 
adopt R v = 3.2 for NGC 6231. This value is 
the same as the total to selective extinction ra- 
tio toward the int ermediate-age open cluster Mil 
(|Sung et alJll999h . 

In addition, the stars near the cluster center 
show somewhat larger color excess ratios, as seen 
in Figure [5] The abscissa is the di stance from the 



cente r of the polarization angle (jFeinstein et al 



2003). But when the abscissa is distance from the 



cluster center, the result does not differ very much. 
The stars at larger distance from the cluster cen- 
ter (r po \ > 14') show very normal ratios, but stars 
in the inner region show large values and a large 
scatter. The Ry of the inner region from 4 optical 
and near-IR bands is 3.43 ± 0.04. SBL98 pro- 
posed that the deficit of low-mass PMS stars with 
Ha emission and the hole seen in reddening map 
are due to the stellar wind from massive stars, WR 
stars or past supernova explosion(s). This higher 
value of the color excess ratios in the cluster center 
may also be related to these effects. 

The extra extinction due to the slightly higher 
value of the color excess ratios of the stars in 
the cluster center is about 0.04 - 0.05 mag in V 
and 0.02 mag in V — I, if we assume the fore- 
ground reddening toward the Sa gittarius arm to b e 
E(B - V) fg = 0.2 - 0.25 mag (ISung et alJl2000h . 
This extra extinction may slightly influence the 
distribution of stars in the middle panel of Figure 
[3] and the upper panel of Figure [7] 

3.3. Color-Magnitude Diagrams and Dis- 
tance 

The color-magnitude diagrams (CMDs) of 
NGC 6231 are presented in Figure [71 From the 
CMDs at the top left we can easily recognize the 
well-developed sequence of early-type MS/post- 
MS members in the left of each CMD. The zero- 
age main sequence (ZAMS) with mean reddening 
and adopted distance is over-pl otted. In addition , 
a few PMS evolution tracks of lSiess et all d200Ctl 
and the 1 M evolution model of iBaraffe et al 



(|1998I ) are superimposed. The empirica l color- 
tempe rature relations for low-mass stars of lBessell 
(Il995l) and synthetic color - temperature relations 
of lBessell et al.l (|1998l ) were employed in the trans- 
formation of physical parameters to observational 



colors. 

Most of the O and early B type stars in NGC 
6231 are X-ray emitters. In addition, many faint 
X-ray emission stars are about 3 mag brighter 
than the adopted ZAMS line, and have the typ- 
ical brightness of PMS stars in young open clus- 
ters. There are not as many X-ray emission stars 
at V — I 1.0, but stars in the lower MS and tran- 
sition zone are mostly X-ray emission stars. The 
stars at V — I ~ 1.0 are field stars at a similar dis- 
tance to NGC 6231. They are relatively old disk 
stars in the Sagittarius arm, and therefore they are 
mostly inactive in X-rays. In addition, the small 
number of X-ray emission stars implies that the 
PMS stars evolve quickly at this stage. 

There are five early-type Ha emission stars with 
V — I < 1.0. Two of them are within the cluster 
radius (see section 3.4 for the radius of NGC 6231), 
the brighter three (ID 294 = Wray 15-1546; ID 357 
= HD 326324 = Hen 3-1269; ID 6017 = ALS 3833 
= Hen 3-1281) are outside the radius of NGC 6231. 
There are two X-ray emission stars at V — I « 0.6 
and V « 10.5 (ID 1499 = SBL 19 = HD 326343, 
B1V; ID 4840 = CD-41 11062B). These two stars 
are more reddened (E(B - V) = 0.67 and 0.71, 
respectively) than the other early-type stars. ID 
4840 is the fainter component of CD-41 11062, and 
is an optical double with ID 4842. 

The distance modulus of NGC 6231 can be 
determined from the distribution of the distance 
moduli of individual early-type stars after cor- 
recting for interstellar reddening. The extra ex- 
tinction mentioned above does not seriously af- 
fect the estimate of distance to NGC 6231 be- 
cause the distance modulus of individual early- 
type stars is estimated using the ZAMS relation in 
the (V, B — V) diagram. The distribution shows a 
peak at V^ — My = 10.6 - 11.0 mag, and only a few 
stars with a larger value. We adopt the distance to 
NGC 6231 as V - M v = 11.0 mag. This value is 
in excellent agree ment with the value obtained by 



Sana et alj (|2005l ). Although there is a slight dif- 



ference in the adopted reddening law, this value 
is identical to the valu e obta ined by SBL98 and 
adopted bv ISana et al. ( 20061 ). The adopted dis- 
tance modulus and mean reddening were applied 
in Figure [3 

We also checked the completeness of our pho- 
tometry in the AV (= V — Vbright) versus distance 
from the nearest bright star. If the brightness of 



6 



a star is brighter than or comparable to the lo- 
cal brightness around a bright star, the brightness 
of the star can be measured. But if the bright- 
ness of a star is fainter than or comparable to 
the local fluctuation of surface brightness around 
a bright star, the brightness of the star cannot be 
measured. And therefore the existence or nonexis- 
tence of a faint star in the photometric catalog is a 
function of the distance from a bright star as well 
as the difference in brightness between two stars. 
The tracing of the avoidance zone makes it pos- 
sible to estimate the area of the incomplete zone 
due to the presence of bright stars. And thus the 
ratio of the total area and the area of the incom- 
plete zone could represent the completeness of the 
photometry. As the seeing is one of important fac- 
tor affecting the completeness of photometry, the 
completeness from this method is for the case of 
bad seeing of 2000, June, 22. The completeness of 
V = 19 mag stars at the cluster center (r < 2') is 
about 60%, but is higher than 80% at r > 2'. Al- 
though the limiting magnitude of our photometry 
is about V = 20.5 mag, the magnitude of over 80% 
completeness, is about V = 19 mag, except for the 
cluster center (r < 2'). The dashed lines in the V 
versus V — I diagram are the apparent loci of PMS 
stars (X-ray emission stars and Ha emission stars) 
of NGC 6231. From the completeness test results, 
our photometry is about 80% complete for PMS 
stars with mass > 0.8 M Q , except for the cluster 
center (r < 2'). 

There is a relatively well-defined sequence of 
foreground stars in the CMDs. They are about 
5.5 mag brighter than the adopted ZAMS of NGC 
6231. Several X-ray emission stars can be found 
in the sequence. The distance to the foreground 
objects is about 130 pc. As NGC 6231 is located 
at the south-eastern edge of the Upper Centaurus- 
Lupus group of the Sco-Cen OB association, the 
stars are likely members of the Upper Centaurus- 
Lupus moving group. 

3.4. Surface Density Profiles and the Ra- 
dius of NGC 6231 

NGC 6231 is the core cluster of the Sco OBI 
association, and therefore the radius of the cluster 
is not easy to determine because many members 
of the Sco OBI association, as well as field stars 
in the Sagittarius arm, are seen in the same sight 
line. In particular, the photometric characteristics 



of Sco OBI association members are very similar 
to those of the cluster members. 

To determine the radius of NGC 6231, we as- 
signed the membership of each star. Stars with 
V < 13 mag and U — B < 0.0 mag were classi- 
fied as massive MS members of NGC 6231. And 
stars in the PMS locus (the upper panel of Figure 
[7|) were classified as PMS members if X-ray emis- 
sion and/or Ha emission is detected, or as PMS 
candidates if no membership information is avail- 
able. Next we determined the center of NGC 6231 
from the geometric center of stars with a certain 
brightness limit. Different brightness limits gave 
a slightly different value, but the difference was 
not large. We calculated the geometrical center 
of the massive MS members with masses greater 
than about 5 M© (7 w 11 mag) and adopted this 
as the center of NGC 6231. The adopted center 
of NGC 6231 is a J20 oo = 16 fc 54 m ll. s 7, <5j 20 oo = 
—41° 50' 20", which is approximately e quidistance 
from HP 152248 C07.5III+0 7III(f) - ISana et all 
(teOOlh l. HD 152249 (09Ib(f) - ISana et all (12008^ . 



and HD 152270 (WC7+06V - Ifflll et a 1.1 (l2()()(lli i. 



We calculated the surface density profiles for 4 
member groups - massive MS/post-MS members 
(V < 10, equivalent ly m > 10 M©), intermediate- 
mass MS stars (V = 10 - 13, equivalently m = 
10 - 2.5 M©), bright PMS stars (V < 15 in PMS 
locus), and faint PMS stars (V = 15 - 18 in PMS 
locus) , and these are shown in Figure [5] In the 
calculation of surface density, we used an annu- 
lus width of V for three faint groups, but a 2' 
width was used for the bright massive group be- 
cause of the small number of massive stars. In 
addition, we also calculated the surface density of 
presumably non-member stars (bright field stars 
- stars brighter and redder than the PMS locus; 
faint field stars - stars fainter and bluer than the 
PMS locus). The error in the surface density is as- 
sumed to follow Poisson statistics, i.e. e = \fN / S 
(where N and S are the number of stars and the 
area of the annulus, respectively). As expected, 
the surface density profiles of field stars did not 
show any radial variation. The results are sum- 
marized in Table |4] 

The fifth column of Table 0] represents the con- 
trast of surface densities between cluster and field. 
The contrast is very high for MS and evolved mas- 
sive stars, but those of PMS stars are very low. 
The surface density of MS/post-MS stars shows 
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a large fluctuation at r > 10' due to the small 
number of early-type stars in the field. Although 
there are large differences in surface density vari- 
ation among member groups, the surface density 
of each group reaches the level of the surface den- 
sity of the field region at r ps 10' regardless of 
the stellar mass. From the radial variation of sur- 
face density profiles, the radius of NGC 6231 was 
adopted to be 10'. 

We also attempted to fit the cluster surface den- 
sity profiles to a Ki ng empirica l density profile of 
a spherical system ( King Tl962T ) . as shown in the 
left panel of Figure |9] As expected, the tidal radii 
of all groups were very similar to each other, and 
between 10' - 20'. But the core radius varied sys- 
tematically from 1.'5 to 5. '5. As King profiles are 
strictly speaking not appropriate for open clus- 
ters, we also fitted the cluster surface density using 
lElson et all (|l987t )'s profiles (EFF profiles). The 
best fit EFF profile parameters are presented in 
Table [4] Using the best-fit values, we computed 
the core radius of the EFF profile, i.e. the radius 
at which the surface density drops to half its cen- 
ter value. As for the King profile, we observed a 
systematic increase from 1.'9 to 4. '7 (correspond- 
ing to 0.9 to 2.2 pc given V - M v = 11.0) for 
populations with decreasing masses. Both King 
and EFF profile fittings imply a mass segregation 
within NGC 6231 - massive stars are more concen- 
trated in the cluster center, while faint PMS stars 
are more widely distributed. We will come back 
to this aspect in Section 4.5. 

4. THE INITIAL MASS FUNCTION 
AND THE AGE OF NGC 6231 

4.1. Adopted Calibrations and Stellar 
Evolution Models 

In SBL98 and our studies on other young open 
clusters, we used the s pectra l type - tempera- 
ture scale of ICrowtherl (|l997h for O type stars. 
Advances in astronomical instruments, observing 
techniques, and non-LTE model atmospheres have 
made an appreciable change in the temperature 
scale of massive O type stars. In this study, 
the spe ctral type - temperat ure scale of O type 
stars of iMartins et al.l (|2005f ) is adopted. Their 
calibration gives somewhat lower temperatures, 
and hence smaller bolo metric corr e ctions for O 
type stars than that of ICrowtherl (|1997h . For 



intermediate- and low-mass stars no significant 
change in the temperature scale was made, and 
therefore the same calibrations were used as in 



Hur et all (120121) . 

In the 2000s, stellar rotation was introduced 
i n the evolutionary m odels of massive stars 
(|Mevnet fc Maederll2000h . B ut in most studies the 
stellar evolution models by ISchaller et aD (|l992h 
were used because new models with stellar ro- 
tation were limited to only a few high masses. 
Recently, new stellar evolution models with rota- 
tion covering a wide mass range were published 
(jEkstrom et al.ll2012t iBrott et al.ll201ll) The stel- 
lar evolution models by IBrott et al.l ( 2011 ) arc 
limited for the stars with masses 5-60 M Q . 
This mass range is useful for NGC 6231, but 
our long term project for investigating the shape 
of the IM F of y o ung massiv e clust ers (e.g. see 



Hur et all (|2012t ): iLim et al.l (|2012b ) 



requires a 

much wide range of stellar masses. In this study, 
we will mainly use th e stellar evolution models by 



Ekstrom et al.l (|2012l ) that include a moderate ini- 



tial rotation (v 



eg 



0Av n 



thus sa 220fcms~ 1 for 



massive main sequence stars). Although the final 
ZAMS temperatures and luminosities of a 1 M Q 
star are very similar t o each other, t he de tailed 
evolutionary trac ks of Baraffe et al.l ([l998) and 
Siess et al ] (|2000h differ from each other as shown 



in Figure El As the models by ISiess et al.l (|2000l ) 
cover a relatively large mass range, we will use 
the PMS evolution models by ISiess et al.l (|2000l) 
for the age and mass estimate of low-mass PMS 
stars. 

4.2. Hertzsprung-Russell Diagram and 
Age of NGC 6231 

The age of an open cluster can be determined 
by comparing the observations with theoretical 
stellar evolution models. The HRD is the basic 
diagram to diagnose a stellar system. To con- 
struct the HRD of NGC 6231, we should correct 
the interstellar reddening of individual stars. For 
early-type stars the reddening is estimated from 
the color-color diagram and applied individually. 
But for low-mass PMS stars, as there is no direct 
way to estimate the reddening from photometric 
data alone, the reddening of these stars are esti- 
mated from the reddening map presented in Figure 
21 Using the adopted calibrations, we transformed 
the reddening-corrected colors to physical parame- 
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ters, such as effective temperature and bolometric 
magnitude. For massive O and early B type stars, 
the spectral type is the most important indicator 
of the effective temperature and hence bolometric 
correction as well. 

The HRD of NGC 6231 is presented in Figure 
[TOl The theoretical ZAMS and two isochrones of 
ages 4 Myr and 7 Myr are superimposed. These 
isochrones were derived by interp olating the the 



oretical stellar evolution tracks of lEkstrom et al 
(2012) and the PMS evolution tracks o fS iessetal 



20001 ). The age of massive stars in NGC 6231 
is between 4 Myr and 7 Myr. This is some- 
what larger than the age obtained by SBL98 
(2.5 - 4 Myr). Although there was a change in 
the stellar evolution models used in the age es- 
timate, the main cause of the big difference in 
age is the change in the adopted calibrations, es- 
pecially for massive O type stars. For another 
young open cluster Westerlund 1 , the new and old 
stellar evolution models from the Genev a group 
(jEkstrom et alj|2012t ISchaller et alJll992h do not 
give any significant difference in age ( Lim et all 
2012 ). However, the stellar evolution models by 



Brott et al.l (|2011l ) give a slightly younger age (3 



4 Myr) for NGC 6231 (see section 5.2 for details). 

The situation for PMS stars is the opposite - the 
age of low-mass PMS stars are younger than pre- 
sented in SBL98. Many PMS members and candi- 
dates are slightly redder than the 4 Myr isochrone. 
The age distribution of low-mass PMS stars is 
shown in Figure 1111 The histograms show the age 
distribution of all PMS stars in the PMS locus 
and are corrected for the contribution of interlop- 
ers (Sco OBI association members and field stars), 
respectively. As the PMS evolution models do not 
give a cons i stent age for the whole mass ran ge 
(|Sung et al.l l2004h (see also (|Sung et al.l Il997h ). 
the distribution is derived only for those PMS 
stars with masses between 0.8 - 1.0 M . The age 
of low- mass PMS s tars is between 1 Myr and 7 Myr 
([Sana et al.l (|2007l) also found a similar spread of 
age of low-mass PMS stars), and therefore the age 
spread is about 6 M yr, which is the typica l spread 
from CMD analysis (|Sung fc BesselfclOj ). 



4.3. The Initial Mass Function 

The mass of an individual star is estimated from 
the HRD. For massive stars it is not easy, because 
a small difference in age gives a very different po- 



sition in the HRD, and it is therefore impossible 
to use a single isochrone or mass-luminosity rela- 
tion. In addition, the complex evolutionary tracks 
of massive stars make this matter even more diffi- 
cult. We tried another way to estimate the mass 
of massive stars - the mean initial mass of stars 
at a given point in the HRD is calculated from 
a Monte Carlo simulation of a stellar system. A 
few million stars with masses larger than 20 M 
were generated for about 500 Myr at a constant 
rate. The stellar masses were generated according 
to an IMF with Salpeter type slope (T = —1.3). 
In the model calculation, stars were evolved ac- 
cording to the recen t stellar evolution models of 



Ekstrom et al.l (|2012T ) . From the model, we calcu- 



lated the mean mass and age of stars within the 
specified photometric box in the HRD. From this 
database we calculated the mass and age of mas- 
sive evolved stars. 

The mass of intermediate-mass MS stars was 
estimated using the mass-luminosity relation of 
stellar evolution models. The mass of PMS stars 
was estimated by interpo l ating the PMS evolu- 



tion models of Siess et al. (2000). And then the 



number of stars in a logarithmic mass interval of 
A log m — 0.2 was calculated. The number of stars 
in the field region was also calculated in the same 
way. After subtracting the contribution of field 
stars, the cluster mass function was then obtained 
by dividing the projected area and mass interval. 
To minimize the effect of binning, we also calcu- 
lated the mass function by shifting the mass bin by 
log m = 0.1. Actually we only estimated the initial 
mass of individual stars, not the current mass of 
the star. Therefore the mass function is the same 
as the IMF except for the most massive bin for the 
clusters with age > 3.5 Myr. For NGC 6231, the 
mass function obtained here is the IMF of NGC 
6231 for m < 45 M Q . 

The IMF of NGC 6231 i s presented in Figure 
HU The IMFs of NGC 2 264 (ISung fc BesselM HOl) 
and the r\ Carina nebula Hur et all (|2012h are also 
presented. The slope of the IMF of NGC 6231 is 
slightly shallower than that of the r\ Carina nebula, 
but the difference is marginal because the shal- 
low slope is partly caused by the dearth of stars 
at logm « 0.2. But the difference between NGC 
6231 and NGC 2264 is evident. The slope of NGC 
2264 is very steep. Although the surface density 
at logm « is very similar, the surface density at 
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logm w 1 of NGC 6231 is about 1 times higher 
than t hat of NGC 2264. Recently iDamiani etaT 



(|2009h obtained a slightly steep IMF (r = -1.3 
- -1.5) for NGC 6231. Such a difference in the 
IMF slope may be caused by the difference in the 
adopted distance, the amount of reddening correc- 
tion, and/or the adopted stellar /PMS evolution 
models. 

4.4. The Total Mass of NGC 6231 



The total mass o f a cluster is an impo rtant pa- 
rameter. Recently IWeidner et all (j2010l ) studied 
the relation between the most massive star and 
its parental cluster mass. They classified clusters 
into three categories according to their total clus- 
ter mass - lowest mass clusters (M c i < 10 2 M ), 
moderate mass clusters (10 2 M < M c \ < 10 3 
M ), and rich clusters (M cl > 10 3 M©). And 
they found a plateau of a constant maximal star 
mass (m max w 25 M ) for clusters with masses 
between 10 3 M Q and 4 x 10 3 M . 

We calculated the total mass of NGC 6231 by 
simulating model clusters using a Monte Carlo 
method. For the model clusters we assumed the 
IMF of NGC 6231 in four cases for the low-mass 
regime - (1) the IMF of NGC 62 31 and that of 
NGC 2264 (ISung fc Besselll l2010h were adjusted 
between logm = 0.6 - 0.4, and we used the ad- 
justed IMF of NGC 2264 for logm < 0.4, (2) the 
IMF of NGC 6231 and that of NGC 2264 were 
adjusted between logm = 0.2 - 0.0, and we used 
the adjusted IMF of NGC 2264 for logm < 0.0, 
(3) the IMF of NGC 2264 was used for logm < 
0.0 without an y adjustment, and (4) the IMF of 
iKroupal (|2002f) . For the massive part of the first 
three cases, we simply extrapolated the IMF of 
NGC 6231 up to logm = 2.0. We assumed the 
upper mass limit of stellar mass as 100 M , and 
an age spread of about 3 Myr (age: 4-7 Myr) 
for the massive stars (m > 30 M ) as obtained 
in section 4.2. A total of 100,000 stars were gen- 
erated, and a scaling factor was calculated to re- 
produce the number of stars with masses between 
logm = 0.4 - 1.0. The total masses estimated 
were 3.0 x 1O 3 M , 2.2 x 10 3 M , 2.1 x 10 3 M , 
and 3.3 x 10 3 M , respectively. And therefore 
the total mass of NGC 6231 is 2.6 (± 0.6) xlO 3 
M . This value is slightly smaller than the total 
mass of NGC 6231 (4595; 2 ^g M ) estimated by 



the total mass may be related to the subtraction 
of field contamination, especially for stars in the 
Sco OBI association. 

4.5. MASS SEGREGATION 

Mass segregation in young open clusters is a 
hot i ssue in the dynamical evolution of star clus- 
ters. iMcMillan et al.l (|2007j ) proposed a plausible 



scenario for mass segregation as mentioned in the 
introduction. NGC 6231 is an important target 
for such a study as the cluster is little obscured, 
has many members, and relatively small differen- 
tial reddening across the cluster. To see the extent 
of mass segregation in NGC 6231, we calculated 
the radial variation of the IMF as shown in Fig- 
ure [13J The IMF for each ring (width of the ring 
= 2') was calculated in the same way as in sec- 
tion 4.3. But as the number of massive stars was 
very small, only one mass bin was used for the 
stars more massive than 10 M . For the outer- 
most ring (r = 8' ~ 10'), the number of stars in 
the ring, after subtracting the field contribution, 
was negative in many cases, and therefore we did 
not try to calculate the IMF for the ring. The 
completeness of photometry was better than 90% 
even in the innnermost ring for logm > 0.4, so 
we did not correct for the effect of data incom- 
pleteness in the left panel of Figure [T3] But the 
effect of incompleteness in logm < 0.4 cannot be 
neglected, especially for the innermost ring where 
the crowding effect is most severe. In the right 
panel of Figure [131 we also calculated and pre- 
sented the slopes of the IMF for a given ring, with 
and without correction for data incompleteness. 

As can be seen in Figure Q21 the slope of the 
IMF increases systematically from the center to 
the peripheral region. The slope of the IMF in 
each ring was calculated, and shows a variation 
with radius from the cluster center (right panel). 
The slope of the IMF varies systematically - it is 
very shallow in the center and very steep at the 
peripheral regions. The right panel of Figure 1131 
shows that the mass segregation in NGC 6231 is 
evident and very systematic. 



Allison et al. 



Weidner et al.l (|2010l ). Such a large discrepancy in 



(2009a) introduced a parameter, 
the mass segregation ratio Amsr, to detect and 
quantify the level of mass segregation in a clus- 
ter. The Amsr compares the path length from 
the minimum spanning tree of a certain kind of 
objects and that of random samples. The merit 
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of this method is that the Amsr is insensitive to 
uncertai nty in determining the center of an open 
cluster. Sana et all (|2010l ) applied this method 
to the young compact cluster in the r\ Carina 
nebula, Trumpler 14, and found a signature of 
mass segregation for bright members (m > 10 
M ), and marginal mass segregation for less bright 
members (K s < 10.5 mag). We also calculated 
the mass segregation ratios for member stars (MS 
stars and PMS stars in the PMS locus within the 
adopted cluster radius) and present it in Figure 
ll"4l In the calculation of Amsr, we used the same 
bin size (N = 20). The most massive two bins 
(m > 8 M Q ) show a distinct mass segregation, 
and stars with masses between 10-3 M Q arc 
marginally mass segregated. It is worth noting 
that as the completeness of photometry is lower 
for the fainter stars, especially fainter stars in the 
center, the Amsr is affected by the data incom- 
pleteness. Given the limitation, we can say that 
massive stars (m > 8 M©) in NGC 6231 are mass 
segregated. Discussion on the origin of mass seg- 
regation is presented in section 5.3. 

5. DISCUSSION 



5.1. HD 153919 

The runaway 06.5Iaf star HD 153919 has 
been proposed to have originated from Sco OBI 
(|Bessell et al.lll974t lAnkav et al.ll200ll ). If the lo- 
cation of HD 153919 U a = 2.80±0 59mas /yr, p, 5 = 
4.71 ± 0.32mas/yr, Ivan Leeuwenl (|2007t) ) before 
the supernova explosion is assumed to be at the 
center of the distribu tion of polarization vectors 
( Feinstein et al. 20031 ). we can check whether or 
not that is a reasonable assumption, on the basis 
of stellar evolution theory. As the angular dis- 
tance between the the polarization center and HD 
153919 is (Aa = l.°92, AS = 3.°84, p = 4.°29), 
the proper motion of the star should have the 
minimum value (/x — e^) for the star to have 
crossed NGC 6231 in the past. And therefore 
if we take the proper motion of the star as 
p! a = 2.21mas/yr, p! & = 4.39mas/yr, then the 
star would have been ejected from the polariza- 
tion center 3.14 ± 0.01 Myr ago. As the age of 
massive stars in NGC 6231 is 4 - 7 Myr, the past 
supernova exploded when the age of NGC 6231 
was 1-4 Myr. 

The life-time of the most massive star (120 M Q ) 



in lEkstrom et all (|2012l) is 3.5 Myr, and that of a 
85 M Q is 4.0 Myr. If HD 153919 originated from 
NGC 6231, the initial mass of the supernova pro- 
genitor would be 85 M Q , or larger. This result is 
consistent with the current stellar evolution mod- 
els and the IMF of NGC 6231, as we could expect 2 
- 3 stars more massive than the most massive star 
in the current NGC 6231. While the end prod- 
uct of the evolution of star typically more mas- 
sive than 25 M Q is expected to be a black hole, a 
massive star may end its life as a neutron stars, 
i f it loses enough mass through wind mass-loss 
(jFrver et al.ll2002l ) and/or binary interaction with 
a lower mass companion. As most massive stars 
are in close binary systems (jSana et al ]l2008lh part 
of the primary's mass is likely to be transferred to 
the companion during the MS stage or in later 
evolutionary stages, and finally it becomes a neu- 
tron star after a supernova explosion. Given that 
HD153919 is a high-mass X-ray binary, the second 
scenario is the more likely. 

5.2. Age of Massive Stars 

The current state-of-the-art stellar evolution 
models do include the effect of stellar rotation. It 
is beyond the scope of this discussion to provide 
a comparison between the different assumptions 
as well as the implementation of different physi- 
cal mechanisms; however, we tested the impact of 
adopting the different available models on the de- 
rived properties of the high-mass stars, and more 
specifically on their derived age. 

Figure [15] (left panel) compares the isochrones 
obtained from the non-rotating evolutionary mod- 
els of Brott et al. and Ekstrom et al. The main 
difference resides in the position of the end of the 
main sequence in the HRD and results from dif- 
ferent assumptions on the value of the core over- 
shooting parameter. The age derived for the mas- 
sive stars is, however, very similar, and a sin- 
gle isochrone, corresponding to an age of about 
3.5 ± 0.5 Myr, reproduces most of the single and 
binary O star population in NGC 6231. 

The age obtained from evolutionary tracks that 
incorporate initial stellar rotation is, however, 
somewhat diffe r ent in both sets of models. The 
Ekstrom et al.l ( 2012 ) rotating models adopt an 



initial v eq of 0.4u cr i t , which is of the order of 200 - 
220 km s _1 for the mass range considered. Figure 
[TBI (right panel) compares the Ekstrom et al. and 
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Brott et al. models with similar initial rotation 
rates. For the Brott et al. models, the location 
of the isochrones is little affected and our previous 
age estimates stand. However, the Ekstrom et al. 
rotating isochrones are shifted to higher tempera- 
tures so older isochrones, of 4 to 7 Myr, are needed 
to reproduce the HRD position of the NGC 6231 
O-type stars. One has to adopt an initial rotation 
larger than 400 km/s for the Brott et al. models 
to fit a similar cluster age. Such a large initial 
rotation is unlikely, however, given that the ma- 
jority of unevol ved O-type stars are relatively slow 
rotators (e.g., (|Pennv fc Giesll2009t )). 

The question of whether the formation of high- 
mass stars in a star forming region triggers or sup- 
presses the fo rmation of lower mass stars (see e.g. 
Herbid (Il962h ) can in principle be tested by com- 
paring the age of the high- and low-mass stellar 
population in NGC 6231. Unfortunately, uncer- 
tainties in the evolutionary calibrations makes it 
difficult to age-date the high-mass stars with suf- 
ficient confidence. 

If the bulk of the NGC 6231 high mass stars are 
truly 5 or 6 Myr old, they would be born at the 
very beginning of the star forming event that lead 
to the formation of NGC 6231. On the other hand, 
an age of 3.5 ± 0.5 Myr old would see them formed 
when star formation activity was at its peak (see 
Fig. QUI. 

5.3. Primordial Mass Segregation Or Dy- 
namical Evolution 

The current theoretical point of view on mass 
segregation in young open clusters is that it 
results from the rapid dynamic al evolution of 
clumps in the subvirial state ([McMillan et al. 
20071 : 1 Allison et al. Il2009bh . This picture implies 



that there should be several clumps in a star form- 
ing region, and that the star forming region may 
be, on average, less dense than the case of only a 
single star forming clump in a star forming region. 
This means that massive stars could be formed in 
the peripheral region of a star forming region. In 
this case, we could have a chance to observe the 
formation process of massive stars with masses 
over 30 Mq, which would be extremely impor- 
tant, as we still we have no detailed information 
on the formation processes of massive stars. 

The surface density profiles and the radial vari- 



ation of the IMF show clearly that NGC 6231 is 
mass-segregated. And the age of NGC 6231 is 
much older than the time required for mass seg- 
regation in a subvirial cluster with substructure 
(see Fig. 2 of lAllison et al. I (|2009bl )). So, is the 
mass segregation in NGC 6231 of dynamical ori- 
gin? Other types of data can help decide. The 
binar y fraction of mas sive O type stars in NGC 
6231 (|Sana et alJl2008h is 100% at the center, and 
it reaches about 60 % from 2' up to 15' (Sco OBI 
region) . The mass of a binary system is more mas- 
sive than single stars, but the binary fraction of O 
type stars is nearly the same, even in the periph- 
eral region of NGC 6231. This fact implies that the 
mass segregation in NGC 6231 may be of primor- 
dial origin rather than the dynamical evolution 



scenari o of iMcMillan et al.l (|2007f) : lAllison et al. 
(|2009bh . 



6. SUMMARY 

New wide-field photometry for the young open 
cluster NGC 6231 is presented and analyzed. The 
main results obtained from this study are summa- 
rized as follows. 

(1) The reddening law toward NGC 6231 was 
derived from color excess ratios from the optical to 
mid-IR Spitzer IRAC bands. The reddening law 
was slightly higher (Ry = 3.2), but very close to 
the normal reddening law. Interestingly, stars in 
the cluster center show relatively larger color ex- 
cess ratios. This fact, and the relatively smaller 
reddening in the center, are conjectured as the re- 
sult of strong stellar winds from massive stars or 
a supernova that exploded a few Myr ago. 

(2) The IMF and age of NGC 6231 were derived. 
The age of NGC 6231 is 4 - 7 Myr (Ar MS = 3 Myr) 
for massive stars, and 1-7 Myr (Atpms = 6 Myr) 
for low-mass PMS stars usin g a new stellar evolu- 
tion models with rotation ( Ekstrom et al.l 1201 2h 
and by adopting new s pectral type - tempe rature 
relation of O type stars ( Martins et al.|[2005 ). The 
slope of the IMF is slightly shallow, but very close 
to the Salpeter value (r = -1.1 ±0.1). In ad- 
dition, the total mass of NGC 6231 is estimated 
to be about 2.6(±0.6) x 10 3 M© from simulated 
model clusters using a Monte Carlo method. 

(3) The surface density profiles of four mass 
groups were derived. Although there was a large 
difference in the ratio of central density to field 
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density, the four groups showed similar cluster 
radii. In addition, the radial variation of the mass 
function was derived and showed an evident mass 
segregation in NGC 6231. The mass segregation 
of massive stars (m > 8 M Q ) was also confirmed 
from the higher value of the mass segregation ra- 
tio Amsr from the minimum spanning tree. The 
mass segregation in NGC 6231 may be of primor- 
dial origin and not as a result of dynamical evolu- 
tion, because of the constancy of the binary frac- 
tion of massive O type stars with radius. 

(4) The runaway O type supergiant HD 153919 
could have originated from NGC 6231. Using the 
Hipparcos proper motion data, the past supernova, 
i f it exploded at the center of the polarization 
(Fcinstein et al ]|2003l) . would have exploded about 
3.1 Myr ago, and the initial mass of the supernova 
progenitor is estimated to be 85 M Q or higher. 
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Fig. 1 — Finder chart of NGC 6231 for the stars 
brighter than V — 17. The size of the dots is 
proportional to the brightness of the star. Large 
open circles represent the position of bright stars 
(data from SIMBAD database). The position of 
stars is relative to the adopted center of NGC 6231 
(a( J2000) = I6 h 54 m 11/70, 5(J2000) = -41° 50' 
20. "3). Large squares represent four fields of view 
of the SSO SITe 2048 x 2048 CCD. 
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This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 
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Fig. 2. — Comparison of the photometry. The 
differe nce A is in the sen se of our new photometry 
minus ISung et al.l (|l998l ) for V < 17 mag (V, V - 
/, and B - V) and 7 < 16 mag ([/ - B). Our 
new photometry is slightly brighter in V, slightly 
redder in V—I and B—V, but systematically bluer 
in B — V for red stars. The difference in U — B 
is somewhat larger and non-linear. See text for 
details. 
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Fig. 3. — The color-color diagrams of NGC 6231. 
(Top panel) the U — B versus B — V diagram for 
V < 16.5 mag. (Middle panel) the B — V versus 
V—I diagram for V < 18 mag. (Bottom panel) the 
Ha index versus V — I diagram for V < 18 mag. 
Thin and thick solid lines in the upper two panels 
represent the intrinsic and reddened ( E{ B — V) = 
0.47 mag) MS relations, respectively. The dashed 
line in the middle panel presents the intrinsic rela- 
tion of giant stars, while the thin solid line in the 
lower panel represents the H a index versus V — I 
relation for normal MS stars ( Park fc Sung 2002), 
but shifted by 0.2 mag in V — I. Crosses, aster- 
isks, squares, and triangles denote, respectively, 
X-ray emission stars, X-ray emission stars with 
Ha emission, Ha emission stars, and Ha emission 
candidates. 
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Fig. 4. — The smoothed reddening map of NGC 
6231. The lines represent the iso-reddening con- 
tours smoothed with the scale length of 1'. The 
line type and thickness represents different amount 
of reddening E(B — V) as shown in the figure. The 
dot represents the early-type stars used in redden- 
ing determination. 
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Fig. 5.— The reddening law of NGC 6231 from the 
optical to mid-IR. The reddening for a given color 
with respect to E(B — V) is shown. The optical 
and mid-IR data are from this study, while near-IR 
data are from 2MASS. The dashed and solid lines 
represent the normal reddening law = 3.1) 
and the reddening law for NGC 6231, respectively. 
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Fig. 6. — Radial variation of color excess ra- 
tios. The radial distance in minutes of arc is 
measured from the c enter of polarization angle 
(|Feinstein et al.ll2003l) . Squares and crosses rep- 
resent stars with Ha and X-ray emission, respec- 
tively. The early-type stars in the field region show 
a nearly normal ratio, but those in the cluster show 
a large scatter. 
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Fig. 7. — The color-magnitude diagrams. (Top 
panel) the V — I versus V diagram. The thick 
solid line represents the reddened ZAMS relation 
at the NGC 6231 distance. T he thin sol i d line s 
arc the PMS evolution tracks of SiessetaL (2000), 
while the thin dashed line overlaid on the thick 
w hite solid line i s the 1 M PMS evolution track 
of iBaraffe et all (|l998l) . The dashed lines in the 
redder part of the reddened ZAMS line represent 
the upper and lower limit of the adopted locus 
of PMS stars in NGC 6231. (Middle panel) the 
B—V versus V diagram. (Bottom panel) the U— B 
versus V diagram. The other symbols are the same 
as Fig. H 




Fig. 8. — Surface density profiles. (Upper) Sur- 
face density profiles of MS stars and stars in the 
PMS locus of NGC 6231. The horizontal bars rep- 
resent the average surface density of field region 
(r > 10'). (Lower) Surface density profiles of the 
field stars above and below the PMS locus of NGC 
6231. The widths of the annuli are V for all groups 
except the most massive group (Ar = 2'). 





Fig. 9. — Surface density profile fitt ings. (Left) 
King profile fits. (Right) EFF profile ([Elson et al 
19871 ) fits. Note that the two figures use different 
units. 




Fig. 10. — The Hertzsprung-Russell diagram of 
NGC 6231. Triangles represent variables. The 
other symbols are the same as Fig. |3l 
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Fig. 11. — The age distribution of stars in NGC 
6231. White and hatched histograms represent, 
respectively, the age distribution of all the low- 
mass PMS stars within the cluster radius, and 
those that have been statistically subtracted for 
the contribution of field stars. 



Fig. 13.— Radial variation of the IMF. (Left) 
The IMF at different radii. (Right) The radial 
variation of the slope of the IMF. The large dots, 
the solid line, and the slope gamma represent the 
radial variation for the case of no correction for 
the data incompleteness, while open circles and a 
dashed line are corrected for the data incomplete- 
ness. The slope of the IMF varies systematically 
from cluster center to the outer region. The IMF 
at the outermost region is not shown because of 
its large uncertainty due to the small number of 
cluster stars. 
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Fig. 12.— The IMF of NGC 6231. Filled 
squares are the IMF of NGC 6231, while open 
circles and open diam onds are that of NGC 2264 
(ISung fc Besselll l2010h and the rj Carina nebula 



Hur et al.l 120121 ) . The solid lines represent the 



least square fit to the IMFs. The slope of the IMF 
(r) is calculated for the stars with logm > 0.2. 



Fig. 14. — The relation between the mass segre- 
gation ratio (Amsr) and the stellar mass. The 
number of member stars in each bin is 20. The 
dot represents the average mass of each bin, and 
the horizontal bar represents the range of stellar 
masses for a given bin. The Amsr of the stars with 
m > 3Mq shows mass segregation, while those of 
low-mass stars have about unity (Amsr ~ 1 - no 
mass segregation). 



Fig. 15. — The upper part of the Hertzsprung- 
Russell diagram. (Left) Non-rotating models, 
(Right) rotating models. 
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Table 1 

Standard Stars, Atmospheric Extinction Coefficients, time variation coefficients, and 

photometric zero points 



Date of Obs. 


standard stars 


Filter 


Extinction Coefficients 


Time Variation 


Photometric Zero Points a 






/ 


0.101 ± 0.007 


0.000 


22.285 ± 0.011 




E5-R c , E5-O c , E5-T° 


V 


0.182 ± 0.004 


-0.0014 ± 0.0008 


22.791 ± 0.009 


2000, Juno, 22 


E6-U, E661, E6-Y, E6-W 


B 


0.294 ± 0.007 


-0.0035 ± 0.0016 


22.640 ± 0.011 




E7-W, E737, E7-c, E7-X 


U 


0.533 ± 0.016 


-0.0056 ± 0.0034 


20.213 ± 0.022 






Ha 


0.101 ± 0.017 




18.067 ± 0.024 






/ 


0.070 ± 0.006 


+0.0063 b ± 0.0011 


22.284 ± 0.008 




E5-R c , E5-O c , E5-T° 


V 


0.134 ± 0.005 


-0.0048 b ± 0.0014 


22.760 ± 0.008 


2000, Juno, 25 


E6-U, E661, E6-Y, E6-W 


B 


0.271 ± 0.006 


-0.0041 b ± 0.0014 


22.622 ± 0.009 




E7-W° , E737 c , E7-c c , E7-X° 


U 


0.506 ± 0.012 


-0.0044 b ± 0.0014 


20.193 ± 0.018 






Ha 


0.092 ± 0.010 




18.072 ± 0.003 



a Photomctric zero points at midnight 
b for UT < 13 h 

c obscrvcd twice - once at meridian, the other at air mass ^ 2. 
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Table 2 
Photometric Data & 



to 



ID 


a J2000 


£.12000 




V 


V 


- / 


B 


- V 


U - B 




ev 


e V -i 




£U-B 




N, 


)bs 




2MASS 


M b D/V c 


SBL98 


Sp 


3610 


16 


:54: 


17. 


95 


-41 


:48 


:01. 


.6 


11.063 





.355 





.210 


-0.492 





.001 





002 





.008 





.005 


2 


2 


2 


2 


16541793-4148013 




480 


B2IVn 


3611 


16 


:54: 


17. 


96 


-41: 


:44 


:57. 


.1 


16.460 


1. 


.256 


1 


.060 


0.654 





.000 





002 





.004 





.021 


4 


4 


3 


2 


16541795-4144568 




6777 




3612 


16 


:54: 


17. 


97 


-42: 


00 


:37. 


.3 


17.808 


1. 


.717 


1 


.345 


0.809 





.013 


0. 


027 





.026 





.108 


1 


1 


1 


1 


16541794-4200373 








3613 


16 


:54: 


17. 


98 


-42: 


Of) 


:07 


.4 


12.898 





.843 





.650 


0.129 





.000 


0. 


003 





.001 





.003 


3 


3 


3 


4 










3614 


16 


:54: 


17. 


98 


-41: 


31 


:22 


.8 


16.849 


1. 


311 


1 


.049 


0.496 





.019 





020 





.022 





.036 


2 


2 


1 


1 


16541799-4131229 








3615 


16 


:54: 


17. 


99 


-41: 


56 


:12 


.5 


17.182 


1. 


.592 


1 


.250 


0.651 





.001 





002 





.017 





.060 


2 


2 


1 


1 


16541798-4156125 




6779 




3616 


16 


:54: 


18. 


01 


-41: 


32 


:47. 


3 


17.420 


2 


.991 


2 


.393 







.005 





014 





.054 






2 


2 


1 





16541801-4132474 








3617 


16 


:54: 


18. 


01 


-41 


:59 


:12 


.7 


17.577 


1 


.667 


1 


.334 


0.873 





.022 





.029 





.027 





.081 


2 


2 


1 


1 


16541798-4159126 




6781 




3618 


16 


:54: 


18. 


03 


-41: 


.45 


:11. 


.8 


14.154 





.907 





.670 


0.381 





.007 





.006 





.005 





.005 


4 


4 


4 


3 


16541802-4145117 




481 




3619 


16 


:54: 


18. 


17 


-41: 


.48 


:13. 


.0 


17.542 


1. 


.414 


1 


.273 







.000 





.011 





.030 






2 


2 


1 





16541819-4148125 




6790 




3620 


16 


:54: 


18. 


17 


-41: 


50 


:16 


.9 


13.065 





.572 





.370 


-0.040 





.002 





.007 





.012 





.000 


3 


3 


3 


3 


16541814-4150164 




482 


B9V 


3621 


16 


:54: 


18. 


18 


-41: 


:53 


:36 


.9 


13.337 





.608 





.420 


0.105 





.011 





.012 





Oil 





.003 


2 


2 


2 


2 


16541819-4153372 




483 




3622 


16 


:54: 


18. 


19 


-41: 


.45 


:33 


.1 


12.094 


2 


.350 


1 


.865 


2.328 





.001 





.017 





.006 





.011 


3 


2 


3 


2 


16541820-4145330 




484 




3623 


16 


:54: 


18. 


22 


-41: 


33 


:59 


.2 


16.877 


1. 


.218 





.971 


0.470 





.005 





Oil 





.009 





.045 


2 


2 


2 


1 


16541823-4133592 








3624 


16 


:54: 


18. 


27 


-41: 


:34 


:12 


.0 


13.276 





.568 





.384 


0.229 





.003 





.005 





.005 





.008 


2 


2 


2 


2 


16541828-4134120 








3625 


16 


:54: 


18. 


29 


-41: 


50 


:10 


.1 


17.168 


1. 


.827 











.061 





.052 










2 


2 








16541824-4150097 




6792 




3626 


16 


:54: 


18. 


29 


-41: 


:51 


:35. 


.2 


9.597 





.360 





.181 


-0.684 





.007 





.012 





.009 





.007 


1 


1 


1 


1 


16541832-4151357 


V 


486 


BlV(n) 


3627 


16 


:54: 


18. 


30 


-41: 


:46 


:38 


.3 


14.640 


1. 


.855 


1 


.629 


1.696 





.015 





.002 





.003 





.019 


4 


4 


4 


2 


16541828-4146381 




485 




3628 


16 


:54: 


18. 


33 


-41: 


32 


:23 


.7 


8.454 





.234 





.116 


-0.456 





.016 





.038 





.024 





.018 


1 


1 


1 


1 


16541837-4132240 






B9Ib-II 


3629 


16 


:54: 


18. 


34 


-41: 


.45 


:00 


.0 


17.447 


1. 


.646 


1 


.258 


0.711 





.022 





.001 





.005 





.147 


3 


3 


2 


2 


16541832-4144599 




6800 




3630 


16 


:54: 


18. 


37 


-41: 


:35 


:23 


.0 


17.462 


1. 


.440 


1 


.175 


0.918 





.008 





.050 





.013 





.082 


2 


2 


1 


1 


16541836-4135228 









a Tablc[2]is presented in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding its form and content. Units of right ascens 
seconds of time, and units of declination arc degrees, arcminutes, and arcseconds. 

b mcmbcrship - X: X-ray emission star, H: Ha emission star, h: Ha emission candidate, + : X + H, -: X + h 

c duplicity or variability - D: stars whose PSF shows a double, but measures as a single star, V: variable from the literature, variable types if available noted in the last column (dSct: c 
(5 Ccphci type variable) 



Table 3 

Comparison with Photoelectric and CCD Photometry 



Author 



Bol^e^d 
Poinstoin 



Scggcwiss (1968 
Schild et al. (1969 
^l^j^son^^^chUd 
Hcskc &c Wcndkcr 
Perry et al. (1990 
Sung et al. (1998 
Piatti et al. (1998 
Baume et al. (1999 




Ay 


n(n cx ) b 


A(V-I) 


n(n cx ) b 


A(B - V) 


n(n ox ) b 


&{U - B) 


n(n ox ) b 


-0.022 ± 0.053 


10 (2) 






+0.009 ± 0.023 


12 (0) 


-0.027 ± 0.045 


11 (1) 


+0.025 ± 0.132 


7 (1) 






+0.017 ± 0.041 


7 (1) 


-0.054 ± 0.061 


7 (1) 


-0.017 ± 0.051 


24 (1) 






+0.007 ± 0.068 


23 (1) 


+0.012 ± 0.091 


23 (1) 


-0.024 ± 0.017 


16 (1) 






-0.003 ± 0.014 


16 (1) 


-0.003 ± 0.030 


17 (0) 


-0.007 ± 0.017 


16 (1) 






+0.000 ± 0.018 


23 (0) 


+0.026 ± 0.040 


23 (0) 


+0.010 ± 0.069 


9 (0) 






-0.017 ± 0.040 


9 (0) 


-0.005 ± 0.040 


8 (1) 


+0.017 ± 0.036 


24 (4) 






+0.013 ± 0.015 


26 (1) 


-0.009 ± 0.035 


26 (1) 


-0.016 ± 0.021 


1157(148) 


+0.027 ± 0.023 


1227(74) 


+0.012 ± 0.018 


1114(173) 


+0.033 ± 0.065 


998(122) 


+0.059 ± 0.056 


58 (9) 


+0.078 ± 0.064 


61 (6) 










+0.010 ± 0.048 


466(102) 


+0.056 ± 0.050 


493(50) 


+0.000 ± 0.055 


496(62) 


+0.018 ± 0.055 


271(60) 



^Comparison made for V < 16 for U ~ B, and V < 17 for V, V - I, and B — V . 
b The number of stars in parenthesis is excluded in the comparison 



Table 4 

Summary of Surface Density Analysis 



Group 


V range 


Pent" 


P/d" 


(Pent - Pfd)/Pfd 


N ste r 




massive MS/post-MS 


< 10 


1.11 


1.34 x 10~ 2 


82.1 


45 




MS stars 


10 - 13 


2.86 


5.00 x 10~ 2 


56.2 


145 




bright PMS 


< 15 


2.23 


3.94 x 10 _1 


4.7 


567 




faint PMS 


15 - 18 


4.46 


1.41 


3.2 


1817 




bright non-mcmbcr 


above PMS locus 




1.68 




1901 




faint non-member 


below PMS locus 




1.74 




1979 




EFF profile fitting parameters 




Po [pc 


a [pc] 


7 


r [pc] 


n [pc d ] 


Ntot 


massive MS/post-MS 


6.80 ± 0.38 


1.63 ± 0.04 


5.45 ± 0.07 


0.88 ± 0.02 


0.048+ 0.005 


32.88 


MS stars 


10.18 ± 1.73 


1.39 ± 0.18 


2.97 ± 0.22 


1.07 ± 0.14 


0.090 ± 0.042 


127.15 


bright PMS 


10.16 ± 1.46 


8.96 ± 1.21 


42.34 ± 0.52 


1.63 ± 0.22 


1.791 ± 0.094 


126.93 


faint PMS 


10.21 ± 1.63 


7.48 ± 3.42 


17.25 ± 1.65 


2.16 ± 0.99 


6.494 ± 0.227 


235.32 



^Surface density at center and field region (r > 10'), respectively. Uint: # / □ 



22 



